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Defects in crystalline silicon consisting of a silicon self-interstitial atom and one, two, three, or four hydrogen 
atoms are studied within density-functional theory (DFT). We search for low-energy defects by starting from an 
ensemble of structures in which the atomic positions in the defect region have been randomized. We then relax 
each structure to a minimum in the energy. We find a new defect consisting of a self-interstitial and one hydrogen 
atom (denoted by {/,H}) which has a higher symmetry and a lower energy than previously reported structures. 
We recover the {/,H2} defect found in previous studies and confirm that it is the most stable such defect. Our 
best {/,H3 } defect has a slightly different structure and lower energy than the one previously reported, and our 
lowest energy {/,H4} defect is different to those of previous studies. 

PACS numbers: 61.05.-a, 61.72.jj, 71.15.Dx, 71.15.Mb 



I. INTRODUCTION 

Hydrogen is a very common impurity in semiconduc- 
tors whose roles in silicon include passivating surfaces and 
defects. 1 Much is known about the vacancy in silicon but 
rather little understanding of self-interstitials has been gleaned 
from experiments, so there is considerable reliance on theoret- 
ical work. Self-interstitials are common in silicon and they are 
expected to react with impurities to form defect complexes. 
Mobile hydrogen atoms are expected to bind strongly to self- 
interstitial defects in silicon, and hydrogen-silicon complexes 
have been detected in experiments 3 

Silicon self-interstitials are readily formed during device 
manufacture and bombardment with electrons or ions. Ac- 
cording to DFT calculations, the most stable structure is 
the split- (110) defect, with the hexagonal interstitial being 
slightly higher in energy and the tetrahedral interstitial being 
still higher in energy. 3 ' 4,5 The results of two quantum Monte 
Carlo calculations are consistent with these three defects hav- 
ing low energies. 3 ' 6 Mukashev et al. 1 attributed the AA12 
electron paramagnetic resonance center to a self-interstitial 
defect, possibly a single self-interstitial. Calculations by Eber- 
lein et a/. 8 found the doubly-positively-charged single self- 
interstitial to be stable at the tetrahedral site, and that it was 
broadly consistent with the AA12 defect. 

Much of what is known about hydrogen in silicon has been 
learnt from studies of vibrational modes which are accessible 
to infrared absorption experiments and may also be calculated 
within first-principles methods. Only one hydrogen-silicon 
complex has so far been firmly identified in experiments. 
Budde et al 2 identified the silicon self-interstitial with two hy- 
drogen atoms using Fourier transform infrared (FTIR) absorp- 
tion spectroscopy. The observed properties were found to be 
in excellent agreement with the results of DFT calculations. 2 

Throughout this paper we denote the silicon self-interstitial 
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atom by / and the n hydrogen atoms by H„, and the whole de- 
fect is referred to as {/,H„}. A metastable defect is indicated 
by an asterisk. We also use the notation devised by Gharaibeh 
et alX where (n)-(m) ■ ■ ■ means n hydrogen atoms bonded to 
one silicon atom and m hydrogen atoms bound to a neighbor- 
ing silicon atom, etc. For example, the silicon self-interstitial 
bonded to two hydrogen atoms mentioned above is referred to 
as {/, H2} and the arrangement of H atoms shown in Fig. |4]is 
described as (l)-(l). 

The first calculations we are aware of on the {/,H} defect 
were by Deak et al.r^- who reported two possible structures 
which, however, are different from those found subsequently 
by Van de Walle and Neugebauer. 1 1 In each structure found 
by Deak et al. 10 the second nearest neighbour silicon atom to 
the hydrogen has a dangling bond, whereas Van de Walle and 
Neugebauer's lowest energy structure does not contain dan- 
gling bonds. Budde et al£ found {/,H} structures based on 
split-(HO) and split-(lOO) self-interstitials both of Cu, sym- 
metry, with the (110) defect being 0.24eV lower in energy. 
The {/,H} defect based on the split- (1 10) interstitial was sim- 
ilar to that found by Van de Walle and Neugebauer 11 whereas 
the split- (100) is similar to the lowest energy structure of 
Deak et a/.— Gharaibeh et al.™ found a structure similar to 
Van de Walle and Neugebauer's 1 1 lowest energy defect. 

Deak et al. 10 ' 12 also studied the {/,H2} defect. Their low- 
est energy {/,H2} defect is also based on a split- (110) self- 
interstitial with its two dangling bonds terminated by hydro- 
gen atoms. Further evidence in favour of this structure has 
been obtained in a variety of studies i 2 ' 9 ' 11 ' 13 Gharaibeh et alX 
also found this to be the most stable of the {/,H„} family of 
defects. 

Hastings et alr^ studied the {/,H3} defect within Hartree- 
Fock (HF) theory, finding a structure with two hydrogen atoms 
bonded to a silicon and a third bonded to a neighbouring sil- 
icon, i.e., a (2)-(l) configuration of hydrogen atoms. They 
found interstitial silyl (SiH3) to be 0.44 eV higher in energy 
and concluded that it is unlikely to be found in bulk silicon. 
More complete DFT calculations by the same group 9 found a 
{/,H3} defect with a (!)-(!)-( 1) configuration. 
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Hastings et al— studied the {/,H4} defect within HF the- 
ory, finding a ground-state structure similar to the {/,H2} 
defect, but with additional hydrogen atoms bonded to the 
nearest-neighbor and second-nearest-neighbour silicon atom 
to the self-interstitial, resulting in a (1)-(1)-(1)-(1) configu- 
ration. They also found a metastable (2)-(2) configuration 
0.2 eV higher in energy and showed that interstitial silane 
(SiH4) is very unlikely to form. A later DFT study by the 
same group 9 found the ground state of {/,H4} to be a (3)-(l) 
configuration. 

In this paper we present calculations for the defects {/, H, }, 
i = 1 ,4, as found by a "random structure searching" approach 
using first-principles DFT methods. We describe the random 
structure searching scheme in Sec. Ill Al The non-standard 
Brillouin-zone integration scheme we have used is described 
in Sec. Ill Bl and the details of our DFT calculations and some 
convergence tests are reported in Sec. Ill CI Our results are de- 
scribed in Sec. [TV] and our main conclusions are summarized 
and discussed in SecM 



II. COMPUTATIONAL APPROACH 
A. Random Structure Searching 

"Random structure searching" has already proven to be 
a powerful tool for finding structures of solids under high 
pressures ! 14,15 ' 16 ' 17 ' 18 The basic algorithm is very simple: we 
take a population of random structures and relax them. This 
approach is surprisingly successful and its performance for 
large systems can be improved by imposing constraints. The 
constraints we have typically employed in work on high- 
pressure phases are to (z) choose the initial positions of a 
local minimum and randomly displace the atoms by a small 
amount, (ii) insert "chemical units" (for example molecules) 
at random rather than atoms, (Hi) search within structures of a 
particular symmetry, and (iv) constrain the initial positions of 
some atoms. 

Pickard and Needs 14 showed that "random structure search- 
ing" can be applied to finding structures of point defects. In 
the current work we have searched for silicon self-interstitial 
defect structures using a 32-atom body-centered-cubic unit 
cell of diamond-structure silicon. The initial configurations 
were generated by removing an atom and its four nearest 
neighbours, making a "hole" in the crystal, and placing six 
silicon atoms at random positions within the hole. The ini- 
tial configurations therefore consisted of a region of perfect 
crystal and a defect region in which atoms are positioned 
randomly. This is an example of a constraint of type (z'v) 
mentioned above. Relaxing the members of an ensemble of 
such initial configurations generated the split- (110), tetrahe- 
dral and hexagonal interstitial configurations, which various 
DFT calculations have shown to be lowest in energy i 2 ' 4 i 5 ' 19 ' 20 

We explored four choices of the hole used to generate the 
initial configurations: (a) remove one silicon atom and its four 
nearest neighbours; (b) remove one silicon atom; (c) do not re- 
move any silicon atoms and take the center of the hole to lie at 
the hexagonal site; (d) the same as (c) but with the hole at the 



tetrahedral site. The initial configurations were generated by 
placing the appropriate atoms randomly within cubic boxes of 
sizes ranging from 2 to 6 A centered on the hole. Choices (b), 
(c), and (d) generally led to us finding the structure we believe 
to be most stable within roughly 100 configurations, while the 
larger hole of choice (a) was less successful and sometimes 
failed to find the ground state within 500 configurations. 

B. Brillouin Zone Sampling Scheme 

The importance of performing accurate Brillouin zone inte- 
grations when calculating defect formation energies has been 
emphasised by, among others, Shim et al— DFT calculations 
by Gharaibeh et al. 9 explored the convergence of BZ sampling 
for self-interstitial-hydrogen complexes in silicon. Instead 
of standard Monkhorst-Pack (MP) sampling 22 we have used 
the multi-k-point generalization of the Baldereschi mean- 
value point scheme 23 outlined by Rajagopal et alM& Three 
linearly-independent reciprocal lattice vectors (bi,b2,b3) are 
chosen, which need not be primitive, and a / x mxn grid of 
k-points is defined by 

/bi jb 2 kb 3 , , „ . 

ki jk =—-\ 1 hk B (7,m,n), (1) 

/ m n 

where 

1 = 0,1,...,/— 1; jf = 0,l,...,m-l; k = 0, 1, . . . ,n - 1. (2) 

The standard Baldereschi mean-value point 23 can be written 
in terms of the b, 

k B = Ofjbi + a 2 b 2 + 0!3b3, (3) 

which defines the a,-, and the offset ks (l,m,n) for the multi- 
k-point scheme 24 ' 25 is 

i, n \ aibl , a2b2 , a3b3 

k B (Z,m,/i) = — h 1 . (4) 

/ m n 

This is the natural multi-k-point generalization of 
Baldereschi's scheme as it corresponds to sampling the 
Baldereschi mean-value point of the supercell obtained by 
choosing direct (real-space) lattice vectors (/ai,ma2,«a3), 
where the a,- are related to the b; by the standard dual 
transformation. We will refer to this as the Multi-B scheme. 

In Fig.Q]we report a comparison for bulk silicon of standard 
nxnxn Monkhorst-Pack grids, the same grids but centered 
on k=0 (Multi-r), and the Multi-B scheme. For odd values of 
n the MP and Multi-r grids are the same, although we could of 
course introduce an appropriate shift of the MP grid for odd 
n which gives smooth convergence. The MP scheme gives 
smaller errors than the Multi-r scheme for even n. The Multi- 
B energies converge smoothly with n and give the smallest 
errors for each value of n. 

The cost of a BZ integration is determined not by n but by 
the number of symmetry inequivalent k-points in the grid and 
the most efficient k-point grid generally depends on the sym- 
metry of the structure. Symmetric structures normally have 
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FIG. 1: (Color online) Magnitudes of the energy differences between 
the converged energy (taken from a Multi-B calculation with n = 20) 
and energies obtained with different k-point grids for a 2-atom cell 
of diamond-structure silicon. The energies are in eV per cell and the 
number of points in each grid is rc 3 . The standard Monkhorst-Pack 
grid, labelled MP, is shown in red [dashed line], Multi-r denotes 
grids centered on k=0 shown in green [dotted line], and Multi-B de- 
notes the multi-k-point generalization of the Baldereschi scheme is 
shown in blue [solid line]. 



more inequivalent k-points in the Multi-B grid than in the cor- 
responding MP grid. During the search stage, symmetry is not 
imposed and we use all n 3 grid points, and therefore the Multi- 
B scheme is the most efficient. We have used the Multi-B grid 
in all searches reported here. 



C. Density-functional theory calculations 

Our calculations were performed using the Generalized 
Gradient Approximation (GGA) density functional of Perdew, 
Burke and Ernzerhof (PBE). 26 The plane-wave basis-set code 
CASTEP— was used with its built-in ultrasoft 28 pseudopoten- 
tials which include non-linear core corrections. 29 All of the 
results presented here were obtained with non-spin-polarized 
calculations. Some searches and large supercell calculations 
were repeated allowing spin-polarization, but no significant 
changes were found in the energy differences between struc- 
tures. 

We carried out convergence tests on the formation energy 
of a bond-centered hydrogen atom in silicon from interstitial 
H2 using 256-atom silicon cells. We chose these systems be- 
cause they contain very different bonding arrangements (Hb c 
contains Si-Si and Si-H bonds while the H2 system contains 
Si-Si and H-H bonds). Hence the estimates of convergence of 
the energy difference between these two systems should give a 
reasonable indication of the convergence of the energy differ- 
ences between other defects. The formation energy, E-p, of a 
bond-centered H in silicon from interstitial H2 calculated with 
iV-atom silicon cells is 



where E(X,N) is the energy of the X defect in a A^-atom sili- 
con cell, and E(N) is the bulk energy for N atoms. 

The Fourier transform grid used for wavefunction manipu- 
lation was set to integrate, without aliasing, frequencies twice 
as high as the maximum frequency in the basis set. We 
checked the convergence of £p(Hbc) with respect to the charge 
augmentation grid required for the ultrasoft pseudopotentials, 
finding it to be converged to within 0.005 eV at 2.75 times the 
maximum frequency in the orbital basis set. We found that 
£p(Hbc) was converged to within ±0.02 eV with a basis set of 
£ , PW = 230eV. 

We tested different k-point sampling schemes in the 256- 
atom cell. The values of Ep(H\ }c ) for the n = 3 standard MP 
and Multi-B grids agreed to within 0.002 eV and we therefore 
considered these converged. The value of ^(Hbc) calculated 
for the 256-atom cell with the n = 2 standard MP grid differed 
from the converged result by ^0.04 e V, whereas the error from 
the Multi-B grid was six-times smaller. 



III. CALCULATING THE FORMATION ENERGIES 

The searches were performed using a body-centered-cubic 
supercell of a size to contain 32-atoms of bulk silicon. We 
used a n — 2 Multi-B k-point grid, which we estimate gives 
energy differences between structures converged to within 
0.006 eV. However, the 32-atom cell is too small to give highly 
accurate geometries and formation energies. We therefore em- 
bedded the most promising structures within 256-atom body- 
centered-cubic unit cells and relaxed using a n = 2 Multi-B 
grid until the forces on each atom were less than 0.001 eVA -1 . 

The formation energy of the self-interstitial is defined as 



E F ({I},N)=E({I},N)-^-E(N), 



(6) 



where E ({/} ,N) is the energy of the self-interstitial cell. 

We define the formation energy per hydrogen atom of a sys- 
tem containing a defect with n hydrogen atoms and i silicon 
atoms relative to a system containing i isolated self-interstitial 
defects and n/2 interstitial hydrogen molecules as 



E F ({I h H n },N) = 



E({Ii,H„},N)-E({Ii},N) 
n 

E{N)-E(H 2 ,N) 



(7) 



£ F (Hbc,A0=£(Hbc,A0 



£(H 2 ,A0 E(N) 



(5) 



Note that Eq. © with i = and n = 1 gives Eq. (|5]l for the 
formation energy of bond centered hydrogen. 

IV. RESULTS 

The formation energies of self-interstitial defects as defined 
by Eq. (O are given in Table |T] In a gree ment with numerous 
previous studies (e.g., Refs. I2TI4I Bi l20t and[l9l) we find the 
most stable defect to be the split- {110} interstitial. The hexag- 
onal and tetrahedral interstitials are 0.03 eV and 0.3 eV higher 
in energy, respectively. Hence all three self-interstitial defects 
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are candidates for forming a low energy {/,H„} defect in the 
presence of hydrogen. 



Defect 


£ F ({/}) (eV) 


Split- (110) 


3.66 


Hexagonal 


3.69 


Tetrahedral 


3.96 



TABLE I: Formation energies for self-interstitial defects in 256-atom 
cells of silicon, as defined by Eq. ((6) 



Our search for hydrogen defects in bulk silicon found the 
lowest energy defects known previously, but no new defects. 
We found the bond centered H atom, the H2 molecule with 
the H-H bond pointing along a (100) direction, and the 
metastable defect of Chang and Chadi. 30 The formation en- 
ergies of hydrogen defects in silicon given in Table [TT] show 
the H2 molecule to be the most stable, as found in previ- 
ous calculations 3 1,32,33 , with the Hj defect and bond-centered- 
hydrogen being respectively O.lOeV and 1.04eV per H atom 
higher in energy. 



Defect 


Configuration 


£p P er H atom (eV) 


Degeneracy dj 


H 2 


Molecule 


0.00 


3 




(D-d) 


0.10 




H bc 


Bond Center 


1.04 





TABLE II: Formation energies per H atom and degeneracies per 
atomic site, dj, of hydrogen defects in silicon. The H2 is the low- 
est in energy defect. 

The formation energies of the hydrogen/silicon self- 
interstitial complexes are calculated with reference to a system 
containing the most stable self-interstitial (the split- (110)) 
and the most stable hydrogen defect in pure silicon, the H2 
molecule. The data in Table Hill shows that the {/,H2} de- 
fect has the lowest formation energy, followed by the {/, H2}*. 
These data are shown in pictorial form in Fig. [2] 



Defect 


Configuration 


Ep per H atom (eV) 


Degeneracy dj 


{AH} 


(1) 


-0.53 


4 


{/,H}* 


(1) 


-0.39 




{/,H 2 } 


(D-d) 


-0.69 


12 


{/,H 2 }* 


(2) 


-0.61 




{/,H 3 } 


(2)-(l) 


-0.39 


24 


U,H 3 }* 


(D-(l)-(l) 


-0.34 




{/,H 4 } 


(l)-(D-(D-(l) 


-0.48 


24 



TABLE III: Formation energies per H atom as defined by Eq. Q 
and degeneracies per atomic site, dj, for various hydrogen/silicon 
self-interstitial complexes. The formation energies are represented 
pictorially in Fig. [2] 

Our searches for the {/,H} defect found the previously- 
known structure^ii of C s symmetry, but we also found a new 
lower-energy structure of higher Ct, v symmetry. The C 3v de- 
fect is 0.14eV lower in energy and is based on a hydrogen 
atom bonded to a hexagonal self-interstitial with the Si-H 



Formation energy per H atom (eV) 



H9 molecule in bulk 




FIG. 2: Formation energies per hydrogen atom using Eq. {7} for var- 
ious hydrogen defects in silicon calculated with respect to silicon 
containing a self-interstitial defect and silicon containing hydrogen 
molecules. 




FIG. 3: (Color online) Two projections of the new {/,H} defect of 
C 3v symmetry found in this work. Silicon atoms are shown in yellow 
(gray) and the hydrogen atom is white. This defect is based on a 
hexagonal self-interstitial rather than the split- (110) self-interstitial 
of previous work. 



bond pointing along a (111) direction, rather than the split- 
(1 10) interstitial on which the C s {/,H} defect is based. 

Our searches found the {/,H2} (l)-(l) defect structure re- 
ported previousl y 2 ' 9 ' 1Q i 1 1 1 1 2 i 1 3 which is shown in Fig. [4] and 
whose formation energy we calculated to be -0.69 eV/H. The 
second most stable {/,H2} defect, which has a (2) structure, 
is the {/,H2}* defect found previously by Deak et alJ£ and 



FIG. 4: (Color online) The most stable {/,H2} defect which has C2 
symmetry and was also found in previous work, 2 . 9 . 10 . 1 '] 12 . 13 Silicon 
atoms are shown in yellow (gray) and the hydrogen atoms are white. 
The defect is based on a split- (110) silicon interstitial defect, with 
hydrogen atoms saturating the two dangling bonds. 




FIG. 5: (Color online) The second most stable {/, H2}* defect which 
has C2 symmetry and was also found in previous work™ Sili- 
con atoms are shown in yellow (gray) and the hydrogen atoms are 
white. The defect is based on a buckled bond-centered silicon self- 
interstitial. 



is shown in Fig. [5] This defect is based on a buckled bond- 
centered Si atom with two H atoms saturating its dangling 
bonds. The {/,H2} and {/,H2}* both have C2 symmetry. 
Hastings et ah 13 found {/, H2}* to be 0.40 eV higher in energy 
than {/,H2} within HF theory. Later Gharaibeh et al. 9 found 
it to be 0.05 eV above the ground state. Our calculations gave 
an energy 0.08 eV/H higher than the {/, H 2 } of Fig. 0] 

The lowest energy {/,H3} defect we found is a (2)-(l) con- 
figuration of C\ symmetry, which is shown in Fig. [6] and has 
a formation energy of -0.39 eV/H. This defect has the same 
configuration of hydrogen atoms but a different structure to 



FIG. 6: (Color online) The new {/, H 3 } defect of C\ symmetry found 
in this work. Silicon atoms are shown in yellow (gray) and the hydro- 
gen atoms are white. The defect is based on a deformed split-(HO) 
self-interstitial. 




FIG. 7: (Color online) The new {/, H4} defect of C\ symmetry found 
in this work. Silicon atoms are shown in yellow (gray) and the hy- 
drogen atoms are white. It is based on the new {/,H 3 } (see Fig. |6) 
with a defect similar to the Hj defect of Chadi 34 adjacent to it on the 
anti-bonding site. 



the one found by Hastings et alJ^ using HF theory. They also 
found a metastable {/,H3}* structure of C\ symmetry with a 
(1)-(1)-(1) configuration just 0.1 eV higher in energy. How- 
ever, more recently this group have used DFT methods and 
found the (1)-(1)-(1) configuration to be lower in energy than 
a (2)-(l) configuration. 9 We find both the new (2)-(l) and (1)- 
(l)-(l) configurations mentioned above, with the new one be- 
ing 0.05 eV/H lower in energy than the (1)-(1)-(1). 

The most stable {/,H4} defect we found, shown in Fig. [71 is 
made up from the {/,H3} defect of Fig.|6]with a defect similar 
to the Hj adjacent to it at the anti -bonding-type site as shown 
by Chadi.-^ This {/,H 4 } defect has a (2)-(l)-(l) configura- 
tion and C\ symmetry. Our {/,H4} defect is different from 
the lowest energy one found by Hastings et alJ^. within HF 
theory, which has a (l)-(l )-(!)-(!) configuration. We have not 
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been able to obtain the lowest energy {/,H4} structure from 
Ref.H although from the description given, we know that it is 
different from ours. We are therefore unable to perform a full 
comparison of energies for {/,H.4}. 

V. DISCUSSION 

We have presented first-principles DFT results using ran- 
dom structure searching for hydrogen/silicon complexes in 
silicon. The searches were carried out in 32-atom silicon cells, 
while the final results were obtained with 256-atom cells. We 
used a multi-k-point generalization of the Baldereschi mean- 
value method to perform the BZ sampling, which we demon- 
strated to be superior to the standard MP sampling. 

Formation energies of the defects were calculated with re- 
spect to the lowest energy self-interstitial defect and lowest 
energy hydrogen defect in bulk silicon. We have confirmed 
that the previously described {/,H2} and {/,H2}* defects are 
the most stable. We have, however, found a new {/,H} defect 
which is significantly lower in energy than the one previously 
reported in the literature. Our defect is based on the hexag- 
onal self-interstitial whereas the previously-reported one was 
based on the split- (1 10) self-interstitial. We also found a new, 
lower energy, {/,H3} defect and a new {/,H4} defect. 

The relative abundances of the defects at zero temperature 
can be calculated from the defect energies as a function of the 
ratio of the concentrations of the self-interstitial and hydrogen 
atoms, «//«h- Fig- Eh shows that only four defects can form 
in this model, {/}, {/,H 2 }, {/,H 4 }, and {H 2 }, with the {/} 
defect corresponding to the lowest-energy split- (110) self- 
interstitial. The main features of Fig. |8]are that when rtj <C «h 
{/,H4} defects are formed and the surplus H atoms form {H2} 
defects, and when «/ ^> «h {AH2} defects are formed and the 
surplus Si atoms form {/} defects. 

At finite temperatures, see Fig. [SJ), we consider the defect 
free energies, which should contain contributions from the vi- 
brational free energy and the configurational entropy. We have 
not evaluated the vibrational free energies, which would in- 
volve very costly phonon calculations, but we have calculated 
the configurational contributions. These are expected to be 
significant because we deal with defects containing from one 
atom (split-(HO) self-interstitial) up to five atoms ({/jtLt}). 
In general, defects containing fewer atoms are expected to 
be favoured by the configurational entropy at higher temper- 
atures. The configurational entropy can be written in terms 
of the number of degenerate defect configurations per atomic 
site, di, where i labels the defect. The degeneracy di can be 
evaluated straightforwardly in some cases. For example, the 



{/,H} defect of Fig. [3] has a degeneracy per atomic site of 
four because it is formed from a Si atom at a hexagonal site 
(of which there are two per atomic site) and the Si-H bond can 
point in one of two directions. Calculating the degeneracy for 
more complicated defects is not necessarily straightforward, 
and therefore we have developed a computational scheme to 
evaluate defect degeneracies. The scheme comprises the fol- 
lowing steps: 

1 . Generate a set of structures by applying the symmetry 
operations of the space group of the host crystal to the 
defect structure; 

2. Identify structures from this set which differ only by a 
translation vector of the lattice of the host crystal and 
remove all but one of them; 

3. The defect degeneracy per primitive unit cell is the 
number of structures remaining. 

The defect degeneracies calculated in this fashion are reported 
in Tables [TT] and [ill] Where we present d\ only for the low- 
est energy defect of each type. The defect concentrations are 
then obtained by minimizing the Helmholtz free energy of the 
system for fixed concentrations nj and «h- In our model we 
consider only the lowest energy defects of each type listed in 
Tables UTI and [TTT1 and in a complete model other defects such 
as clusters of self-interstitial Si atoms and complexes involv- 
ing other impurity atoms should be considered. 

It is interesting to note that only the {/,H2}, {/,H2}*, and 
{/jtLt} defects are perfectly saturated, i.e., each Si atom has 
four covalent bonds and each H atom has one, and that these 
defects have the lowest total formation energies (defined as 
nEp({I,H n }) with E-p given by Eq. (0). It is, of course, not 
possible to achieve perfect saturation of a Si/H structure with 
an odd number of H atoms, as such a structure would contain 
two bonds per Si atom and half a bond per H atom. 

Overall we conclude that "random searching" is a useful 
tool for finding the structures of low-energy defects in semi- 
conductors. 
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